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We present a new derivation of the Casimir force between two parallel plane mirrors at zero 
temperature. The two mirrors and the cavity they enclose are treated as quantum optical networks. 
They are in general lossy and characterized by frequency dependent reflection amplitudes. The 
additional fluctuations accompanying losses are deduced from expressions of the optical theorem. A 
general proof is given for the theorem relating the spectral density inside the cavity to the reflection 
amplitudes seen by the inner fields. This density determines the vacuum radiation pressure and, 
therefore, the Casimir force. The force is obtained as an integral over the real frequencies, including 
the contribution of evanescent waves besides that of ordinary waves, and, then, as an integral over 
imaginary frequencies. The demonstration relies only on general properties obeyed by real mirrors 
which also enforce general constraints for the variation of the Casimir force. 



I. INTRODUCTION 

An important prediction of quantum theory is the exis- 
tence of irreducible fluctuations of electromagnetic fields 
in vacuum. Besides their numerous observable conse- 
quences in microscopic physics, vacuum fluctuations also 
have observable effects in macroscopic physics, for exam- 
ple the Casimir force they exert on mirrors [l|]. 

Casimir calculated this force in a geometrical configu- 
ration where two plane mirrors are placed a distance L 
apart and parallel to each other, the area A of the mirrors 
being much larger than the squared distance A^S> L 2 . He 
considered the ideal case of perfectly reflecting mirrors 
and obtained an expression which, remarkably, depends 
only on the geometrical quantities A and L and on the 
fundamental constants h and c 

Hctt 2 A 

Fcas = 2mA (1) 

This attractive force has been observed in a number 
of 'historical' experiments 0, 0, 0, 0JII_ which confirmed 
its existence and main properties 0, H 0- Several re- 
cent experiments reached an accuracy in the % range 
by measuring the force between a plane and a sphere 
fwi fTTL \l2\ U| or two cylinders 0]. Similar experi- 
ments were also performed with MEMS [THl ITfij | (see also 
\vf\). An experiment studied the plane-plane configu- 
ration considered by Casimir but, as a consequence 
of the difficulties associated with this geometry, reached 
only a 15% accuracy (see reviews of recent experiments 
in HHH). 

The Casimir force is the most accessible experimental 
consequence of vacuum fluctuations in the macroscopic 
world while vacuum energy is known to raise a serious 
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problem with respect to gravity and cosmology (see ref- 
erences in [^1E2)- This is a reason for testing the predic- 
tions of Quantum Field Theory concerning the Casimir 
effect with the greatest care and accuracy. The theory of 
the Casimir force is also a key point for the experiments 
searching for the new weak forces predicted by theoretical 
unification models to arise at distances between nanome- 
ter and millimeter [21E1E1E1E1E1E1,E3 The 
Casimir force is indeed the dominant effect between two 
neutral objects at /im or sub/im distances so that an ac- 
curate knowledge of its theoretical expectation is as cru- 
cial as the precision of measurements in such experiments 

m 

In this context, it is essential to account for the dif- 
ferences between the ideal case considered by Casimir 
and the real experimental situation. Recent experiments 
use metallic mirrors which show perfect reflection only 
at frequencies below their plasma frequency. They are 
performed at room temperature, with the effect of ther- 
mal fluctuations superimposed to that of vacuum fluc- 
tuations. In the most accurate experiments, the force 
is measured between a plane and a sphere, and not be- 
tween two parallel planes. The surface state of the plates, 
in particular their roughness, should also affect the force. 
A large number of works have been devoted to the study 
of these effects and we refer the reader to E3] f° r a 
bibliography. 

The evaluation of the Casimir force between imperfect 
lossy mirrors at non zero temperature has given rise to a 
burst of controversial results [H El El El IM El El 
EM E3> EH which constitutes a part of the motivations 
for the present work. For the sake of comparing experi- 
mental measurements and theoretical expectations, it is 
necessary to have at one's disposal a reliable expression 
of the Casimir force in the experimental situation. In 
the present paper, we focus our attention on the effect 
of imperfect reflection of the mirrors. Other effects, in 
particular the effect of temperature, will be addressed in 
follow-on papers. 

We consider the original Casimir geometry with two 
perfectly plane and parallel mirrors. Except for these as- 
sumptions, we consider arbitrary frequency dependences 
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for the mirrors which, in particular, may be lossy. We 
evaluate the Casimir force as the effect of vacuum radi- 
ation pressure on the Fabry-Perot cavity formed by the 
two mirrors. The net force results from the balance be- 
tween the repulsive and attractive contributions associ- 
ated respectively with resonant or antircsonant frequen- 
cies. It is obtained as an integral over the axis of real fre- 
quencies, including the contribution of evanescent waves 
besides that of ordinary waves. It is then transformed 
into an integral over imaginary frequencies by using phys- 
ical properties fulfilled by all real mirrors. 

The formula obtained here for the Casimir force turns 
out to be identical to the expression already published in 
|42fl but the new derivation has a wider scope of valid- 
ity than the previous one since it remains valid for lossy 
mirrors. The fact that the formula keeps the same form 
despite the widening of the assumptions is intimately re- 
lated to a theorem which relates the spectral density of 
the fields inside the cavity to the reflection amplitudes 
seen by the same fields. This theorem was demonstrated 
in 0| and ^3 m specific cases and we prove it in the 
present paper without any restriction. To this aim, we 
introduce a systematic treatment of lossy mirrors and 
cavities as dissipative networks [J}. We define scatter- 
ing and transfer matrices for elementary networks like 
the interface between two media or the propagation over 
a given length in a medium. We then deduce the matri- 
ces associated with composed networks, like the optical 
slab or the multilayer mirror. 

The results obtained in this manner are therefore ap- 
plicable to a large variety of mirrors, still with the as- 
sumption of perfect plane geometry. In the particular 
case of a slab with a large width, the Lifshitz expression 
|45l l4(| is recovered. At the limit of perfectly reflectors, 
the ideal Casimir formula Q is obtained. More generally, 
the expression gives the Casimir force as an integral writ- 
ten in terms of the reflection amplitudes characterizing 
the two mirrors. This integral is finite as soon as the am- 
plitudes obey the general properties of scattering theory 
already alluded to. In other words, the difficulties usu- 
ally associated with the infinitencss of vacuum energy are 
solved by using the properties of real mirrors themselves 
rather than through an additional formal regularization 
technique. 

We finally show that the same physical properties con- 
strain the variation of the Casimir force. In particular, 
they invalidate proposals which has been done for 'taylor- 
ing' the force at will by usin g m irrors with specially de- 
signed scattering amplitudes |47l l48| . In these proposals, 
the balance between attractive and repulsive contribu- 
tions to the force is change, leading to the hope that the 
Casimir force could reach large or have its sign changed 
from an attractive force to a repulsive one |47| . Using the 
simple model of a one-dimensional space, it has already 
been shown that these hopes cannot be met for ar- 
bitrary mirrors built up with dielectric layers. Here, the 
argument is generalized to the Casimir geometry in three- 
dimensional space with the following conclusions : the 



Casimir force cannot exceed the value obtained for per- 
fect mirrors, it remains attractive for any cavity length 
and its value is a decreasing function of the cavity length. 
This is true for any mirror obtained by piling up layers of 
media described by dielectric functions. This definition 
of multilayer dielectric mirrors includes the case of metal- 
lic layers, provided that magnetic effects play a negligible 
role in the optical response. 

II. VACUUM FIELD MODES 

As explained in the Introduction, we consider in this 
paper the original Casimir geometry with perfectly plane 
and parallel mirrors aligned along the directions x and 
y. This configuration obeys a symmetry with respect 
to time translation as well as transverse space trans- 
lations along these directions. We use bold letters for 
two-dimensional vectors along these directions and de- 
note r = (x, y) the transverse position. As a consequence 
of this symmetry, the frequency lo, the transverse vector 
k = (k x ,k v ) and the polarization p = TE,TM are pre- 
served throughout the scattering processes on a mirror 
or a cavity. The scattering couples only the free vac- 
uum modes which have the same values for the preserved 
quantum numbers and differ by the sign of the longitu- 
dinal component k z of the wavevector. 

In the present section, we introduce notations for the 
vacuum field modes, first in empty space and then in 
a dielectric medium. These notations are chosen to be 
well-adapted to the symmetry of the problem. 

A. Vacuum modes in empty space 

In empty space, the components of the wavevector are 
given for each field mode by the frequency u>, the inci- 
dence angle 9 and the azimuthal angle ip 

k x — |k|cos^) |k| = — sin# 

c 

fcj, = |k|sin<^ k z = — cos 6* (2) 

|k| is the modulus of the transverse wavevector and the 
longitudinal component k z may be expressed in terms of 
the preserved quantities lj and k 



^=^v^~ k2 0=±1 (3) 

4> is defined as the sign of cos 9 and represents the di- 
rection of propagation with +1 and —1 corresponding 
respectively to rightward and leftward propagation. 

The two polarizations p = TE, TM are defined by the 
transversality with the incidence plane of electric and 
magnetic fields respectively. They are given by the unit 
electric vectors e~ 

e x = cos 9 cos ip e x = — sin ip 
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~TM 



cos a sin <p 



e™ = -sin< 



TE 
V 

TE 



= cosy 

= 



or, equivalently, the unit magnetic vectors ft™ = e 
and /3 TE = — e™. For each mode, the wavevector and 
polarization vectors form an orthogonal spatial basis. We 
have chosen linear polarizations described by real compo- 
nents; hence the unit vectors <?and (3 are not affected by 
the complex conjugation appearing below in the relation 
between positive and negative frequencies. 

The two modes corresponding to the same values of ui, 
k and p but opposite values of <f> are coupled by scattering 
on a mirror. For this reason, we introduce a label m = 



(4) 



PTE 



(uj, k, p) gathering the values of u>, k and p. A mode freely 
propagating in vacuum is thus labeled by m and <j) and 
the summation over modes is described by the symbols 
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Note that <fi appears implicitly as the sign of k z in the 
first form whereas it appears explicitly in the second one. 
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The free vacuum fields are then written as linear superpositions of modes 

'hid 



E(v,z,t) = ^/cZ vac 

m<f> 



B(r,z,t) = 



Z. 



VIC 

rrwf) 



2 

hw - 



i(tjt— k.r— k z 



i(a»t— k.r— k z z) 



— i(ujt~\i..r—k z z) i / <fr \ t i(ujt — k.r— k z z) 



(6) 



The vacuum impedance Z vac = Hoc ~ 377f2 describes the electromagnetic constants in vacuum. In the following, the 
symbol e will be reserved to the relative permittivity with the value 1 in vacuum. 



The quantum field amplitudes and (e^j) corre- 
spond to positive and negative frequency components. 
They fit the definition of annihilation and creation op- 
erators of quantum field theory and obey the canonical 
commutation relations 50] 



(27r) 3 <^ (k-k')S(k z -k' z ) 8 PP ,8 4 
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In the vacuum state, the anticommutators of quantum 
amplitudes are derived from the corresponding commu- 
tators 



„e, 



t 



= 



(8) 



The dot symbol represents a symmetrized product. 



B. Stress tensor in empty space 

The energy density per unit volume Too is a quadratic 
form of the fields E and B 



Too (r,z,t) 



1 

2cZ v 



(E 



2 B 2 ) 



(9) 



When subsituting the expression of free fields, Too is ob- 
tained as a bilinear form of the field amplitudes. Here, 
we study the averaged radiation pressure in the vacuum 
state which leads to a contraction m! — m in the sums 
over modes. Using the vacuum property ijSJl. we find the 
averaged energy density in vacuum equal to the sum over 



the modes of 



huj 



(T 00 (r,M)>, 



E 



hid 
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(10) 



As it is well-known, this energy density is infinite. 

The radiation pressure on plane mirrors oriented along 
xy directions is determined by the component T zz of the 
Maxwell stress tensor 



T zz (r,z,t) 



1 



2Z V 



(E-E + c 2 B- B) 



(11) 



Here, the dot symbol represents a symmetrized product 
of the quantum amplitudes and, simultaneously, a scalar 
product of the vectors; the overline symbol describes the 
mathematical reflexion of a vector with respect to the 
plane xy 



E x = E x 



E u 



En, 



E z 



-E z 



(12) 



As for Too, averaging T zz in vacuum state leads to a 
contraction over the modes with the result 

(T Z z ( r j z , £))vac = Z\2 "T" (^m - ^™ 0m-P™J 

711 ip 
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TljJ 



E-7T cos2 ' 
2 

rruj) 



(13) 



This expression is similar to the expression (I10|) of the 
energy density with an extra factor cos 2 9 well-known in 
studies of radiation pressure. The sum over modes is still 
infinite but this infiniteness problem will be solved in the 
forthcoming calculation of the Casimir force. 



C. Fields in dielectric media 

In the following, we consider mirrors built up as dielec- 
tric multilayers. Each dielectric medium is characterized 
by a relative permittivity e [u>] or, equivalently, an in- 
dex of refraction n [to] — y/e [to] depending on frequency. 
The magnetic permeability is kept equal to its vacuum 
value since this corresponds to all experimental situa- 
tions studied so far. We stress again that this definition 
of dielectric mirrors includes the case of metals as long 
as the magnetic response plays a negligible role. We con- 
sider layers thick enough so that the dielectric response 
is local, i.e. described by a wavevector-independent per- 
mittivity e [to] . 

We will sometimes take the plasma model as a first 
description of metallic optical response 



e [u 



LO 2 



LUp 



2-kc 



T- (14) 



where top and Ap represent respectively the plasma fre- 
quency and the plasma wavelength. This simple model 
is not sufficient for an accurate evaluation of the Casimir 
force between real mirrors |3lj |. To this aim, it is nec- 
essary to describe the optical response of metals with 
a dissipative part associated with electronic relaxation 
processes. As a consequence of causality, the real and 
imaginary parts of n are related to each other through 
the Kramers-Kronig dispersion relations [5l| . 

For any function of frequency more generally, causal- 
ity is unambiguously characterized in terms of analyticity 
properties : n [u>] or e [u>] are analytical functions of u) in 
the 'physical domain' of the complex frequency plane, 
that is the domain of frequencies uj with a positive imag- 
inary part 3w > 0. This property is obeyed by other re- 
sponse functions to be encountered below and it will play 
an important role in the derivation of the Casimir force. 
We will introduce an equivalent notation £ for complex 
frequencies with the physical domain now defined by a 
positive real part for £ 



w = i£ 5i£ > 



(15) 



The dispersion relation is changed inside a refrac- 
tive medium to 



k x = |k| cos ip |k| = n [u>] — sin# 

c 

k z = n [u] — cos 9 (16) 
c 



The preservation of u> and k at the traversal of an inter- 
face is equivalent to the Snell-Descartes law of refraction. 

The sign has to be carefully chosen when extracting the 
square root to express k z in terms of the conserved quan- 
tities lo and k. As soon as the refractive index contains 
an imaginary part, this is also the case for k z and the de- 
phasing exp(ik z z) associated with propagation includes 
an extinction factor. In order to ensure that this factor 
is effectively a decreasing exponential, we have to choose 
a specific root defined differently for the two propagation 
directions 6 = ±1 of the field 



k z — i(j)K 



em 



k 2 



5Rk > 



(17) 



k y = |k| sin tp 



The argument has been presented for freely propagating 
modes but it holds as well for evanescent waves confined 
to the vicinity of an interface between two media. In this 
case, the sign of k z is also chosen so that it corresponds to 
an extinction when the distance to the interface increases 
and this choice is still described by equation i|17|) ■ In the 
following, we will use systematically the notations £ and 
k, keeping in mind that the causality relations have to 
be written for each value of the conserved quantity k. 

Besides the dispersion relation l|l(j|) . the dielectric 
medium also changes the impedance, that is the ratio be- 
tween magnetic and electric field amplitudes. Precisely, 
the impedance is changed from the value Z wac in empty 
space to the value — in a dielectric medium of index 
n, resulting in reflection at the interface. 



III. MIRRORS AS OPTICAL NETWORKS 

We now introduce the description of mirrors as optical 
networks. We present the scattering and transfer repre- 
sentations and the relations between them. The transfer 
approach is well adapted to the composition of networks 
which are piled up. We first consider elementary net- 
works such as an interface or propagation inside a refrac- 
tive medium. We then use the composition law to study 
composed networks such as the slab and multilayer. In 
the present section, we only consider classical fields or, 
equivalently, mean quantum fields. The next section will 
be devoted to the full quantum treatment including the 
addition of noise associated with the losses inside the 
mirror. 



A. Scattering and transfer representations 

We first introduce the scattering and transfer rep- 
resentations for an arbitrary network represented with 
two ports and four fields. These fields are identi- 
fied as lefthand/righthand (symbols 'L' and 'R'), right- 
ward/leftward (arrows — > and <— ) or input/output fields 
(labels 'in' and 'out'), as shown on Figure ^ 
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> 

S, 



out < 



out > 



in < 



FIG. 1: Scattering and transfer representations of a network. 



Let us emphasize that the arrows are a symbolic repre- 
sentation of the two modes coupled by the network which 
correspond to the same label m and to the two opposite 
signs <f) = ±1. The geometrical directions of propaga- 
tion are given by the wavevectors of equation (T5J. The 
coupling between the fields is described by reflection and 
transmission amplitudes represented below by scattering 
or transfer matrices. 

In the scattering point of view, we gather the input and 
output fields in twofold columns related by a S— matrix 



\e° 



£\ 



= S\£ 



cout 
|cout\ _ / <-% 
\ L I I £out 



s = 



r t 
t r 



(18) 



r and r are the reflection amplitudes while t and t are the 
transmission amplitudes. We will also use an equivalent 
convention where the output ket is defined with the upper 
and lower components exchanged 



|£ out ) 



cout 

C R 

cout 



S\£ in ) 



v |£ out ) 

S = i]S 



7] = 

t r 
r 1 



1 

1 



(19) 



This convention simplifies some algebraic manipulations 
while being completely equivalent to the former conven- 
tion. For comparison with previous works, note that the 
former notation (|18f) was used in whereas the latter 
one (|19|l was used in |49| . 

In the transfer point of view, the network is de- 
scribed by lefthand and righthand columns related by 
a T— matrix 



\£l) = T \Sb 



\£* 



T = 



£i 



a b 
c d 



(20) 



The matrix rj introduced in (|19fl exchanges the two 
directions of propagation. We also use in the following 
the matrices tt± which project onto each direction 



7T+ 



1 





1 



(21) 



These matrices obey simple rules which define an alge- 
braic calculus in the space M.2 (C) of 2 x 2 matrices with 
complex coefficients 




7T_|_ 7T_ = 7T_ 7T_|_7r_ = 7T_ 7T_|_ = 

rj A = I rjTT + = 7T_7? 777T_ = 7T + ?7 (22) 

The identification of Figure is written as 

IT- \£ R ) = 7T_ |£ in ) 

n.\£ L )=TT-\£^) (23) 

It relates the transfer and scattering amplitudes. We de- 
compose the scattering equations (|18f) on the two com- 
ponents and use (|23[) to rewrite them as 

7T+ |£r) = K + S (?T + \£ L ) + 7T_ \£ R )) 

tt_ \£ h ) = ir_S(ir+ \£ h )+ir_ \£ R )) (24) 
This linear system may be put under a matrix form 



7T_ - SV+ ) \£ L ) 



7T + — SlT- ) \£ 



(25) 



It is equivalent to the transfer equation l|2(J|l with the 
T— matrix obtained as 



T = — ( 7T_ — Sttj 



ir + — Sir- 



(26) 



The converse transformation is obtained by perform- 
ing the same manipulations in the reverse order. Starting 
from the transfer equation (|2U|I and using J22J), one ob- 
tains a linear system which is equivalent to the scattering 
equation 119|) with 



S = — (7T_ — Tn + ) {n + — T7T_) 



(27) 



The relations (|26|l and (|27|l have the same form. They 
represent an idempotent homographic transformation 
in the space M.2 (C), care being taken for the non- 
commutativity of multiplications in this space. When in- 
verting algebraically the homographic relations (|26|l and 
(|27H . one obtains equivalent expressions 



S=(7T + +7T_T) (7r_+7T + iy 

T = (ir + +TT-S) [ir-+ir + S 



(28) 



Other equivalent expressions are obtained from the equal- 
ities 



^7r_ — Sir+j (7r_ — Tir + ) = I 

7T_ + TT + S ) (7T_ + 7T + T) = / 



(29) 



All these expressions may be written in terms of the 
scattering and transfer amplitudes 
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d = 



if 



t = 



t 

ad — be 



(30) 



The more formal homographic transformations written 
above are nevertheless useful, as it will become clear in 
forthcoming calculations. 



B. Composition of optical networks 

The T— matrices are perfectly adapted to the compo- 
sition of optical networks corresponding to a piling up 
process (see Figure |5J|. 



A 


B 











AB 



FIG. 2: Composition of networks : two networks labelled A 
and B are piled up to build up a network AB. 



On each network, the transfer equations are written as 

\£ L {A}) = T{A}\£ K {A}) 

\£ h {B}) = T{B} |£r{B}) (31) 

The brackets { } specify the network for which the 
T— matrix or field column is written. Identifying the 
fields according to Figure © 



|£ L {AB}) ee \£ h {A}) 
\£ R {AB}) ee |£r{B}} 



\£ R {A}) ee \£ L {B}) 



(32) 



we deduce that the piling up process is equivalent to the 
product of T-matrices 



|£ L {AB}) = T{AB} |£ R {AB}) 
T{AB} = T{A}T{B} 



(33) 



We have assumed the two networks to be in the imme- 
diate vicinity of each other but without any electronic 
exchange between them, which again corresponds to the 
assumption of thick enough layers. 



C. Elementary networks 

We now study two elementary networks, that is the 
traversal of an interface and the propagation over a given 
length inside a dielectric medium. 



For the scattering at the plane interface between two 
media with indices no and ri\ , we write the reflection and 
transmission amplitudes as the Fresnel scattering ampli- 
tudes Reflection amplitudes r p {Int} are obtained 
from characteristic impedances z v defined for plane waves 
with polarization p in each medium and from the conti- 
nuity equations at the interface 



r p {Int} 



TE 



TM 



-r p {Int} 

Tli COS 9\ 

no cos #o 

Til COS #o 



1 



z 



p 



1 + Z'P 

Hi 

eiK 



no cos 9\ £oKi 
Then the transmission amplitudes are obtained as 



(34) 



/^{Int} = J-i P {Int} 

Ko V Ki 



1 - (rP{Int}) 2 (35) 



We deduce the expression of the transfer matrix 

1 



T p {Int} = 



e 2 



K V 2 sinh ft? 



e 2 
e 2 



f3 p = In 



z p + l 
zp-1 



(36) 



We now consider the process of field propagation over 
a propagation length £ inside a dielectric medium charac- 
terized by a permittivity e. For this elementary network, 
the T— matrix has the simple form 



T{Prop} = 



k 2 I 



(37) 



The optical depth a does not depend on the polarization. 

Note that the composition is commutative within the 
class of interfaces or that of propagations : it corresponds 
to the multiplication of the z— parameters for interfaces 
and to the addition of a— parameters for propagations. 
But the composition is no longer commutative when in- 
terfaces and propagations are piled up. 



D. Reciprocity theorem 

We now prove a reciprocity theorem obeyed by arbi- 
trary dielectric multilayers, i.e. networks obtained by 
piling up interfaces and propagations. 

To this aim, we first remark that the ratio of the two 
transmission amplitudes is related to the determinant of 
the T— matrix 



= ad — be = det T 



(38) 



This follows from the relations l|3l)|l between S— and 
T— amplitudes for an arbitrary network. Then, it is clear 
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from (|33|) that the determinant of T is simply multiplied 
under composition 

detT{AB} = detT{A} detT{B} (39) 

For the two kinds of elementary networks studied previ- 
ously (see eas l36l37|) . the determinant of T is the ratio of 
the values of K at the right and left sides of the network 

detT=— (40) 

It follows that this relation is valid for any optical net- 
work composed by piling up interfaces and propagations. 

In the particular case where the network has its two 
ports corresponding to vacuum, which is the case for a 
mirror, the values of k are equal on its two sides and the 
T— matrix has a unit determinant 

detT=l t = t (41) 

Note that reciprocity corresponds to a symmetrical 
S— matrix and has to be distinguished from the spatial 
symmetry of the network with respect to its mediane 
plane which entails r = r. 

This theorem is the specific form, when the symmetry 
of plane mirrors is assumed, of the general reciprocity 
theorem demonstrated by Casimir |54| as an extension 
to electromagnetism of Onsager's microreversibility the- 
orem [^1 . We have disregarded any static magnetic field 
which could affect these reciprocity relations. 

E. Slabs and multilayers 

We now consider the dielectric slabs and multilayers 
as composed networks and we deduce their transfer and 
scattering amplitudes from the preceding results. 

The slab is obtained by piling up a vacuum/matter 
interface with indices Hq = 1 and n\ at its left and right- 
hand sides, propagation over a length £ inside matter, 
and a matter/vacuum interface with now n\ and n = 1 
at its left and righthand sides. We denote T{Int} the 
T— matrix associated with the first interface and obtain 
the T— matrix associated with the second interface as the 
inverse of T{Int}. As a consequence of the composition 
law (|33|) . the T— matrix associated with the slab is ob- 
tained as 

T{Slab} = Tiln^TiPmpjTilnt}- 1 (42) 

Using the expressions I|3ul37[l of T{Int} and T{Prop}, 
we evaluate T{Slab} as 

T{ Slab} = -±- ( Sinh & + . , ) (43) 

sinh/3 V - Slnna smb. [ft - a) J v ; 

We deduce the form of the S— matrix which is simultane- 
ously reciprocal (t = t) and symmetrical in the exchange 
of its two ports (r = r) 

„ rm n 1 / — sinh a sinh/3 \ , ... 

5{Slab} = sinh(^ + a) ( sinh^ -sinhaj ^ 



In the limiting case of a small thickness a — > 0, we find 
t{Slab} — * 1 and r{Slab} — > 0, which means that the 
slab tends to become transparent. In this case indeed, 
the propagation can be forgotten and the two inverse 
interfaces have their effects cancelled by each other. 

The opposite limiting case of a large thickness is of- 
ten considered since it fits the usual experimental situ- 
ations. More precisely, experiments are performed with 
metallic mirrors having a thickness much larger than the 
plasma wavelength. This is why the limit of a total ex- 
tinction of the field through the medium is assumed in 
most calculations. This corresponds to the so-called 'bulk 
limit' with e~ Q — * and r{Slab} — * — = r{Int} in 
eg . (|44|l : the reflection amplitude is determined entirely 
by the first interface. Let us emphasize however that 
the bulk limit raises several delicate problems. First, the 
transmission amplitude i{Slab} vanishes in this limit so 
that the T— matrix is not defined, with the drawback of 
invalidating the general method used in the present pa- 
per. Then, the bulk limit cannot be met in the case of 
non absorbing media where e~ a remains a complex num- 
ber with unit modulus for any value of I. Even in the 
presence of absorption, a large value of the width i does 
not necessarily imply a large value of the optical thick- 
ness a since k may go to zero at normal incidence and 
zero frequency, leading to a transparent slab in contrast 
with the results of the bulk limit. Therefore a reliable 
calculation must consider the experimental situation of 
mirrors with a large but finite thickness. In the present 
paper, we consider the general case of arbitrary mirrors 
and test the reliability of the bulk limit in the end of the 
calculations. 

We can deal with the case of dielectric multilayers sim- 
ilarly. If we consider as an example the multilayer ob- 
tained by piling up a vacuum/matter interface with in- 
dices no = 1 and n\ at its left and righthand sides, prop- 
agation over a length l\ inside the medium 1, an interface 
between media 1 and 2, propagation over a length I2 in- 
side the medium 2, and an interface between medium 2 
and vacuum, its T— matrix is obtained as the product 

T{Multilayer} = T{Int01}T{Propl}T{Intl2} 

xT{Prop2}T{Int20} (45) 

Alternatively, the same multilayer may be obtained by 
piling up two slabs each corresponding to one of the layers 

^{Multilayer} = T{Slab010}T{Slab020} (46) 

In the last two equations, the indices specify the different 
interfaces, propagations or slabs using an obvious conven- 
tion. 

Since any multilayer mirror is obtained by piling up 
slabs connecting two vacuum ports and thus obeying the 
reciprocity relation t = t, we can use a simple form of 
the composition law written in terms of scattering am- 
plitudes [49] 

t\ r B _ _ r A tl 
tab = r A + r AB = r B + 

1 - r A r B 1 - r A r B 
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£ab 



1 - r A r B 



(47) 



For readibility, we have specified the networks by using 
subscripts rather than brackets. We will proceed simi- 
larly in forthcoming specific computations. Iterating this 
composition law, we can compute the scattering ampli- 
tudes for any dielectric multilayer. This systematic tech- 
nique is quite similar to the classical computation tech- 
niques used for studying multilayers j5|| . It is generalized 
to the full quantum treatment in the next section. It also 
leads in the following to general results constraining the 
variation of the Casimir force for arbitrary dielectric mir- 
rors. It reproduces the known results for the multilayer 
systems which have already been studied [ljj, |57j • 



IV. QUANTUM TREATMENT OF LOSSY 
MIRRORS 

Up to now, we have performed a classical analysis 
which is not sufficient for the purpose of describing the 
scattering of vacuum fluctuations. Real mirrors consist 
of absorbing media which scatter incident fields to spon- 
taneous emission modes and reciprocally scatter fluctu- 
ations from noise modes to the modes of interest. The 
S— matrix calculated previously cannot be unitary for a 
lossy mirror but it should be the restriction to the modes 
of interest of a larger S— matrix which includes the noise 
modes and obeys unitarity. In the present section, wc 
characterize the additional fluctuations for a lossy mir- 
ror by using the corresponding 'optical theorem', that is 
also the unitarity of the larger S— matrix (see [5^, 
and references therein). 

We assume that the scattering restricted to the modes 
of interest still fulfills the symmetry of plane mirrors 
considered in the previous classical calculations. This 
amounts to neglect multiple scattering processes which 
could couple different modes through their coupling with 
noise modes. Except for this assumption, we consider 
arbitrary dissipative media and discuss the optical the- 
orem in the scattering and transfer points of view. We 
use the latter one to deal with composition of additional 
fluctuations when lossy mirrors are piled up. 



A. Noise in the scattering approach 

Should we use the previous classical equations for the 
quantum amplitudes, we would find that the output fields 
cannot obey the canonical commutators, except in the 
particular case of lossless mirrors. This implies that the 
input / output transformation for quantum field must in- 
clude additional fluctuations superimposed to the classi- 
cal equations 



|e out ) = S*|e in ) + \F) 



(48) 



same matrix as previously and \F) is a twofold column 
matrix describing the additional fluctuations. All these 
quantities depend on the quantum number m which is 
common to all fields coupled in the scattering process. 



FIG. 3: Representation of a dissipative network, with addi- 
tional fluctuations coming from the noise modes. 



The additional fluctuations are linear superpositions of 
all modes coupled to the main modes ef^ by the micro- 
scopic couplings which cause absorption. As an example, 
the atoms constituting a dielectric medium couple the 
main modes to all electromagnetic modes through spon- 
taneous emission processes, represented symbolically by 
the wavy arrows on Figure |3J The stationarity assump- 
tion implies that only modes having the same frequen- 
cies are coupled. In particular, it forbids parametric 
couplings which could couple modes with different fre- 
quencies and 'squeeze' the vacuum fluctuations |6fJ. The 
whole scattering matrix which takes into account all cou- 
pled field modes is unitary and this basic property makes 
the canonical commutation relations compatible for input 
and output fields. In contrast, the reduced scattering ma- 
trix containing only the classical scattering amplitudes 
coupling the main modes is not unitary, except in the 
particular case of lossless mirrors. 

In order to write the unitarity property of the whole 
scattering matrix, it is convenient to represent the ad- 
ditional fluctuations \F) by introducing auxiliary noise 
modes |/) and auxiliary noise amplitudes gathered in a 
noise matrix S' 



\F)=S'\f) 



S' = 



r' t' 
t' f' 



(49) 



O out\ 



\e—) and |e in ) are defined as in (|18|l with the quantum 
amplitudes e in place of the classical fields £, S is the 



The components of the twofold column |/) are defined to 
have the same canonical commutators as the input fields 
in the main modes. In fact, they are linear superpositions 
of the input vacuum modes responsible for the fluctua- 
tion process. They are defined up to an ambiguity : any 
canonical transformation of the noise modes leads to an 
equivalent representation of the additional fluctuations, 
which corresponds to a different form for the noise am- 
plitudes while leading to the same physical results at the 
end of the computations. 

For any of these equivalent representations, the norm 
matrix S' S' ^ has the same expression determined by the 
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optical theorem, that is the unitarity condition for the 
whole scattering process, 

S& + S'S' t = / (50) 

where / is the 2x2 unity matrix. This is easily proven 
by a direct inspection of the explicit expressions of the 
commutators of the output fields. The same inspection 
shows that noise modes corresponding to different values 
of to are not correlated to each other. Condition (|50|l is 
made more explicit when SS' and S' S' ' are developed 
in terms of scattering amplitudes 

rr+rr+rr + 1 1 = tt + rr + 1 1 + r r 

= 1 

rt* + tr* + r't'* + t'r'* = tr* + ft* + tV* + r't'* 

= (51) 

More detailed discussions are presented for the case of 
the slab in appendix 1X1 

The description of noise may as well be represented 
with the alternative representation 1)19(1 of the scattering 
process 

p*) = S\e ia ) + \F) 

j^^j" = 77 |e out > \F)=r)\F) (52) 

The additional fluctuations are then represented in terms 
of the same noise modes and of a modified noise matrix 

\F) = S> \f) S> = rjS' 

S'S' f =I-S& (53) 

B. Noise in the transfer approach 

We now present the description of additional fluctu- 
ations in the transfer approach. Performing the same 
manipulations as in the previous section, we transform 
equation l(52() into 

(tt_ - Stt+) |e L ) =-(*+- Stt_) |e R ) + \F) (54) 

We thus get transfer equations with additional fluctua- 
tions described by a twofold column |G) 

\e L )=T\e R ) + \G) 

|G> = (55) 

The T— matrix has the same expression 1(26(1 as previ- 
ously and the additional fluctuations |G) are a linear ex- 
pression of the fluctuations \F) defined in the scattering 
approach. This linear relation may be written under al- 
ternative forms by using the relations 1(29(1 

\F) = (7r_-^7r + )^> = (^-T7T + )- 1 |G) 

\G) =(n_-Tn + )\F) (56) 



In the scattering approach, the norm of additional fluc- 
tuations is described by matrices S'S' t and S'S' which 
are themselves determined by the optical theorem l(5U|) or 
((53(1 . In order to translate these properties to the trans- 
fer approach, we rewrite 1(55(1 in terms of the canonical 
noise modes |/) and of noise amplitudes gathered in a 
matrix T' 

\G)=T'\f) 

T' = (tt_ - SV+) S' = (tt_ - Ttt+) §' (57) 
The associated norm matrix is 

T'T' t = (tt_ - Tir + ) S'S' 1 (tt_ - Tir + )^ (58) 

Using equations l(53|) and ((27(l . we rewrite it as 

T'T' f = TW* - $ $ = 7T+ - tt_ (59) 

$ is a diagonal matrix with two eigenvalues representing 
the directions of propagation <p = ±1 of the field. 



C. Composition of dissipative networks 

Using these tools, we now write composition laws for 
the fluctuations and their norms. 

We start from transfer equations written for each net- 
work A and B 

|e L {A}) = T{A}|e R {A}) + |G{A}) 

|e L {B}) = T{B}|e R {B}) + |G{B}) (60) 

Using the identifications ((32|) associated with the compo- 
sition law, we deduce for the composed network 

|e L {AB}) = T{AB} |e R {AB}> + |G{AB}> 
|G{AB}> = |G{A})+T{A}|G{B}> (61) 

The fluctuations |G{AB}) are a linear superposition of 
fluctuations |G{A}) and |G{B}) added in A and B. 

In order to obtain the composition law for the norm 
matrices, we develop the additional fluctuations |G{AB}) 
on the canonical noise modes associated with the two 
elements 

|G{AB}) = T'{A} |/{A}) + T{A}T'{B} |/{B}) (62) 

Since the noise modes associated with different elements 
are uncorrelated, |G{AB}) may be rewritten in terms 
of new canonical noise modes and new noise amplitudes 
such that 

|G{AB}) = T'{AB} |/{AB}) 
T'{AB}T'{AB}^ = T'{A}T'{A}^ 

+ T{A}T'{B}T , {B} t r{A} t (63) 
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Using expression (|59|l of the optical theorem for both 
networks A and B, we deduce that the composed network 
AB obeys the same relation 



T'{AB}T'{AB}' ! 



= r{A}$r{A}t - $ 

-T{A} (T{B}$T{B} t - $) T{A}t 
= r{AB}$T{AB}t - $ (64) 



Equivalently, the S— matrix of the composed network AB 
obeys the optical theorem l|50[) as soon as the two net- 
works A and B do. 



D. Resonance for cavity fields 

We have studied the scattering or, equivalently, the 
lefthand/righthand transfer of fields by a composed net- 
work AB. We want now to characterize the properties 
of the fields inside the cavity formed between A and B. 
This problem will play a key role in the evaluation of the 
Casimir force (see next section). 



< 
s, 



> 



< 

Sr 



FIG. 4: Cavity formed within a composed network : L and R 
denote the fields at left and right sides of the network whereas 
C denotes the cavity fields. 



The situation is illustrated by Figure 0] which, in con- 
trast to Figure|21 keeps the trace of the intracavity fields. 
In algebraic terms, the cavity fields are defined by rewrit- 
ing the identifications l-il'l) as 

|e L {AB}> = |e L {A}> |e R {AB}) = |e R {B}) 
|e c {AB}) = |e R {A}) = |e L {B}> (65) 

From now on, we drop the label {AB} for the composed 
network and use subscripts for the networks A and B. 

In order to express the cavity fields in terms of the 
input modes and additional fluctuations, we first write 
the cavity fields |ec) in terms of the righthand ones |e R ) 

|e c ) = T B |e R ) + \Gb) (66) 

We then identify the two components of |e R ) as 



7r+|e R ) = 7r + |e°«t) ^ + (S|e m ) + |F)^ 
7T_ |e R ) = tt_ |e in ) (67) 



Using the expression of \F) in terms of \G) and the com- 
position law (|61|l for |G), we deduce 



i?|e in ) + i?^|/A)+i? B l/B) 



R = T B N 



N 



TT + S + 71 

P = -Nn^ 
R' B = (I + T B PT A ) T B 



(tt_ + n+ry 



R' A = T B PT' A 



(68) 



As already explained, the unitarity of scattering entails 
that the output fields have the same commutators as the 
input ones. But this is not the case for the cavity fields 
which have their commutators determined by the matrix 



Q = RRX 



R' a R' a 



R' B R B 



(69) 



Expanding this quadratic form and using the composition 
law ()63[1 . we rewrite Q as 

Q = TuNN^TX + T B PT'T' ^P^Tl 



+T B PT A T B T^ 



(70) 



Using relation l|59|l for the three networks A, B and AB, 
we obtain a simpler expression after a few rearrangements 



g = -$ - t b pta$ - ^tIp^t^ 



(71) 



We now proceed to explicit calculations of these matri- 
ces. We note that P = —tn + where t is the transmission 
amplitude of the network AB and deduce 



TbPTa* = t 



a B a A 
c B a A 



-a B b A 
-c B b A 



(72) 



t is simply the inverse of the transfer amplitude a asso- 
ciated with the network AB (see ea l30|l and the latter is 
deduced from the composition law (|33f) 



1 

t = - 
a 



a = a A a B + b A c B 



(73) 



Then, the transfer amplitudes of the networks A and B 
may be substituted by the associated scattering ampli- 
tudes, leading to 



T B PT A $ 



1 



1 - ^A^B 



1 r A 
r B r A r B 



(74) 



Collecting these results and proceeding to slight rear- 
rangements, we finally get 



Q = I 



1 



1 - r A r B 



1 



(1 - r A r B ) 



r A r B r A 
r B r A r B 

r A r B r A 
r B r A r B 



(75) 



In the following we will use the diagonal terms of the 
matrix Q to evaluate the Casimir force. 
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E. Scattering on a Fabry-Perot cavity 

In order to prepare the evaluation of the Casimir force, 
we generalize the preceding expression to the case of the 
Fabry-Perot cavity containing a zone of field propagation 
between the two mirrors Ml and M2 (see Figure |5J). 



> 



> 



< 

S, 



< 



> 

s 



<- 



FIG. 5: Representation of a Fabry-Perot cavity : L and R 
denote the fields at left and right sides of the cavity whereas 
C denotes the cavity fields inside the Fabry-Perot cavity; these 
cavity fields are defined at an arbitrary position between the 
two mirrors. 



The distance between the two mirrors is denoted L 
and the cavity fields are defined at an arbitrary position 
inside the cavity, say at distances L\ from Ml and L 2 
from M2 with L\ + L 2 = L. In these conditions, the 
study of the Fabry-Perot cavity is reduced to the problem 
studied in the preceding subsection through the following 
identifications : the network A contains the mirror Ml 
and the propagation LI with Ta — ImiTli while the 
network B contains the propagation L2 and the mirror 
M2 with Tb = Tl,27m2- The transfer amplitudes for the 
networks A and B are derived from those corresponding 
to Ml and M2 and from phase factors corresponding to 
the propagations LI and L2 



t A = t ie ~ ai 
r ie - 2Ql 



a\ = k Lx 
r A = n 

a 2 = k L 2 



= r 2 



r 2 e 



-2a 2 



(76) 



We have labeled the amplitudes for the mirrors Ml and 
M2 with mere indices 1 and 2; kq is defined in vacuum. 
These results entail that the reflection amplitudes ta and 
tb are seen from a point inside the cavity as the product 
of phase factors by the reflection amplitudes r 1 and r 2 
seen from a point in the immediate vicinity of Ml and 
M2. 

We then deduce the scattering amplitudes for the 
whole cavity 



n 



t\r 2 e~ 2a 
D 



r 2 



rit 2 e~ 



t = t = 



t\t 2 e 



D 

D = 1 — r\r 2 e 
and the expression of Q 

g = i + i 

D 



-2a 



ol\ + a 2 



(77) 



1 

~D~* 



r\r 2 e 2a rie 2ai 
r 2 e~ 2a2 T\r 2 e~ 2a 

f\r 2 e~ 2a Tie~ 2ai 
r 2 e~ 2c * 2 r\r 2 e~ 2a 



(78) 



The diagonal terms in the matrix Q coincide with the 
Airy function 



.9=1 + / + /* = 



1 



\r\r 2 e 



|1 - rir 2 e 



/ = 



r\r 2 e 



1 - r x r 2 e 



-2a 



(79) 



D 



This result will play the central role in the derivation of 
the Casimir force in the next section. It means that the 
commutators of the intracavity fields are not the same as 
those of the input or output fields. They correspond to 
a spectral density modified through a multiplication by 
the Airy function g. This is the basic property used in 
Cavity Quantum ElectroDynamics |olT |. 

It is clear from the present derivation that this result 
has a quite general status : it is obtained for any inner 
field in any composed network, assuming the symmetry 
of plane mirrors. This pro perty was already known for 
non absorbing mirrors |42| and for lossy mirrors sym- 
metrical with respect to their mediane plane 01 ■ The 
present derivation proves that it is also valid for arbitrary 
dielectric multilayers with dissipation. The final result 
only depends on the reflection amplitudes r\ and r 2 of 
the mirrors as they are seen from the inner side of the 
cavity. The reflection amplitudes seen from the outer 
side and the transmission amplitudes do not appear in 
expressions (|78I79|) . This can be interpreted as resulting 
from the unitarity of the whole scattering processes. 



CASIMIR FORCE BETWEEN REAL 
MIRRORS 



We may now deal with the radiation pressure of vac- 
uum fields on the mirrors of a Fabry-Perot cavity. We 
show that the resulting Casimir force is a regular integral 
which can be written over real or imaginary frequencies. 
We then derive general constraints obeyed by the Casimir 
force for arbitrary dielectric mirrors. 
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A. Vacuum radiation pressure 



If we first consider a mirror isolated in vacuum, the radiation pressure is obtained by adding the contributions of 
the 4 fields coupled in the scattering process 

( P )vac = J2 ^™ C ° S2 9m ^ e m L ' e m L e rn L ' e m L g rn R ' C R ~ e ™ R ' e m ft) ( 80 ) 
— \ / vac 

m 

The identification of these fields is given by Figure 0). We have developed the sum over </> and kept the symbol m to 
represent the quantum numbers (w, k,p). We assume that the whole system is in vacuum, that is at zero temperature, 
so that the anticommutators of input fields are given by relation |jSJ. Since the commutators are the same for the 
output and input fields, the vacuum radiation pressure vanishes in the case of an isolated mirror. In other words, the 
two sides of the mirror play equivalent roles so that no mean force can appear. 

When we consider two mirrors forming a Fabry-Perot cavity, the two sides of a given mirror are no longer equivalent 
since one is an inner side and the other an outer side. It follows that the compensation observed for an isolated mirror 
does no longer hold, resulting in the appearance of the Casimir force. In order to evaluate the force, we write the 
mean radiation pressures (Pi) vac and (P2) vac on mirrors Ml and M2 (see Figure^ 

( P l)vac =Ys nWm 
m 

( p 2) vac = fiWm 



cos 2 1 



'm ( e m L ' e m L + C m L ' e m L e m C ' e m C e m C ' e m C / 

\ / vac 



COS 2 8 m (emC' e mC +e mC' e mC~ e mR' e mR~ e mR' e mR/ (^1) 



/ vac 



For the same reasons as previously, the field anticom- 
mutators are given by (JSJ for input and output fields. 
For intracavity fields, they are multiplied by the Airy 
function l|79|) like the commutators 



J t 

-m C 



1 r 

2 

1 



C> e m C 



T.9 



m u mm 



(82) 



As shown in the previous section, these expressions do 
not depend on the position inside the cavity where the 
cavity fields are defined. We finally deduce the mean 
radiation pressures on mirrors Ml and M2 



(Pi 



= -uv 



hw m cos 2 m (l- g m ) 



(83) 



At this point, it is worth emphasizing that we have as- 
sumed equilibrium at zero temperature for the whole sys- 
tem : not only the input fields but also any fluctuations 
associated with loss mechanisms inside the mirrors cor- 
respond to zero-point fluctuations, whatever their micro- 
scopic origin may be. Otherwise, the expression of the 
force discussed in the following would be affected. 

The pressures have opposite values on the two mirrors 
Ml and M2. This entails that the global force exerted 
by vacuum upon the cavity vanishes, in consistency with 
the translational invariance of vacuum. In the following, 
we denote F the Casimir force calculated for Ml when 
considering the limit of a large area A 3> L 2 



A (Pi) 



A 



(l-Sm) (84) 



The sign conventions used here are such that the positive 
value obtained below for F corresponds to an attraction 
of the two mirrors to each other. 



B. The force as an integral over real frequencies 

We now perform a change of variable to rewrite the 
summation symbol as specified in JSJ 



F = 



A (Pi 

A J2 



d 2 k 

4^ 



2^ 



hk z (l-s£M) (85) 



We will now specify the domain of integration for u>. 

Up to now, we have discussed the scattering for ordi- 
nary waves which freely propagate in vacuum and corre- 
spond to frequencies w larger than the bound c |k| fixed 
by the norm of the transverse wavevector. But we must 
also take into account the contribution of evanescent 
waves which correspond to frequencies to smaller than 
c |k|. These waves are fed by the additional fluctuations 
coming from the noise lines into the dielectric medium 
and propagating with an incidence angle larger than the 
limit angle. They are thus transformed at the interface 
into evanescent waves decreasing exponentially when the 
distance from the interface increases. As is well known 
|53|. the properties of these evanescent waves are con- 
veniently described through an analytical continuation 
of those of ordinary waves. This analytical continuation 
can only be dealt with in terms of functions having a well 
defined analyticity behaviour. This is not the case for the 
Airy function [u>] but we know that this function is the 
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sum (I79|l of parts having well denned analyticity proper- 
ties 

I 1 - Pk Ml 



/ k p M 



(86) 



p£ [w] is the 'open loop function' corresponding to one 
round trip of the field inside the cavity and defined as the 
product of the reflection amplitudes r£ 1 [w] and r£ 2 [w] 
of the two mirrors and of the propagation phaseshift 
e -2« L. ft j g an ana lytical function in the physical do- 
main of complex frequencies 5i£ > with the branch 
of the square root chosen so that tfin > 0. Since the 
transverse wavevector is spectator throughout the whole 
scattering process, analyticity is defined with k fixed. 

Then, /£ [w] is the 'closed loop function' built up on 
the open loop function p£ [w]. It is also an analytical 
function, thanks to analyticity of the open loop and to 
a stability property which has a natural interpretation : 
the system formed by the Fabry-Perot cavity and the 
vacuum fluctuations is stable because neither the mir- 
rors nor the vacuum would have the ability to sustain an 
oscillation. In some cases, the stability can be derived 
from a more stringent passivity property |49| which may 
essentially be written |p£ Ml < !• However, the passiv- 
ity property is sometimes too stringent to be obeyed by 
real mirrors (see more detailed discussions in appendix 
IbI . In any case, the stability property, i.e. the absence 
of self sustained oscillations, is sufficient for the present 
derivation of the Casimir force. 

We are now able to give more precise specifications of 
the domain of integration in Ij85(l . Using the decomposi- 
tion l|86|l , we write the contribution of ordinary waves to 
this integral as the sum of two conjugated expressions 



ord 



J o: 



d 2 k 

4^ 



|k| 



(111.' 
2^ 



hk z fl M (87) 



The integral J-" or( j is built on the retarded function [oj] 
which may be extended through an analytical continua- 
tion from the sector of ordinary waves to that of evanes- 
cent waves. The contribution of evanescent waves to the 
force is thus obtained as 



J Q. 



(88) 



The final expression of the Casimir force is the sum of 
the contributions of ordinary and evanescent waves that 
is also the integral over the whole axis of real frequencies 



d 2 k 

4^2 



duj . 
2^ 



Kino fl 



(89) 



As far as ordinary waves are concerned, this corresponds 
to the intuitive picture where the Casimir force results 
from the radiation pressure of vacuum fluctuations fil- 
tered by the cavity ji^. The contribution of evanes- 
cent waves is but the extension of the domain of inte- 
gration to the whole real axis with the cavity response 
function fP [lu] extended through an analytical continua- 
tion. In the evanescent sector, the cavity function [lu] 
is written in terms of reflection amplitudes calculated 
for evanescent waves and exponential factors correspond- 
ing to evanescent propagation through the cavity. This 
means that it describes the 'frustration' of total reflec- 
tion on one mirror due to the presence of the other. This 
explains why the radiation pressure of evanescent waves 
is not identical on the two sides of a given mirror and, 
therefore, how evanescent waves have a non null contri- 
bution to the Casimir force. 



C. The force as an integral over imaginary 
frequencies 

Using the Cauchy theorem, we now rewrite the Casimir 
force H89f) as an integral over the axis of imaginary fre- 
quencies. 



Imco=Rei; 



c. 



c. 



c„ 



c|k| 



Reto—Imi; 



F 



F P , 



FIG. 6: Contour representing the frequencies of interest for 
the evaluation of the Casimir force : Co and C c correspond to 
the real frequencies associated with ordinary and evanescent 
waves; Ci correspond to the imaginary frequencies and Coo to 
a quarter circle with a radius allowed to go to infinity. 



Since Ko/k K] l& analytical in the domain 5i£ > 0, its 
integral over a closed contour lying in this domain has to 
vanish. We choose the contour drawn on Figure which 
consists of the positive part of the real axis including 
ordinary (C ) and evanescent (C ) waves, a quarter of 
circle Coo with a very large radius and, finally, the imag- 
inary axis Cj run from infinity to zero. Now the function 
K o/k K] S oes to zero for large values of the frequency, 
as a consequence of transparency at high frequency, a 
property certainly valid for any realistic model of optical 
mirror. Thanks to this property, the contribution to the 
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integral of Coo vanishes. We then deduce that the inte- 
grals over the real axis [0,+oo[ and over the imaginary 
axis [0, +ioo[ are equal. 

We thus get a new expression of the force F as an 
integral over imaginary frequencies u>, that is also as an 
integral over real values of £, 



F 
T 



T 



V 



T* = 2T 
d 2 k 
4^2 



2tt 



tiKo fl [if] 



(90) 



We have used the fact that T is real, so that T* is simply 
equal to T . This property is less obvious, but also true, 
with T written as an integral over real frequencies. We 
wish to emphasize more generally that expression (|90ll is 
mathematically equivalent to (|89|l . The former expres- 
sion is closer to the physical intuition whereas the latter 
is better adapted to explicit computations of the force. 

Expression (|90|) gives the Casimir force between real 
mirrors described by arbitrary frequency dependent re- 
flection amplitudes. It is a regular integral as soon as 
these amplitudes obey the physical assumptions used in 
the derivation : causality, unitarity and high frequency 
transparency for each mirror, stability of the system 
formed by the two mirrors and the scattered vacuum 
fields. The demonstration holds for dissipative mirrors 
and not only for lossless ones. 

The limit of perfect mirrors is obtained in expression 
(|90|l by letting the reflection amplitudes go to unity, 
which leads to the Casimir formula (JIJ . This can be con- 
sidered as an alternative demonstration of the Casimir 
formula without any reference to a renormalization or 
regularization technique. Basically, the properties of real 
mirrors, in particular their high frequency transparency, 
are sufficient to provide a regular expression of the force, 
as it was guessed a long time ago by Casimir Q . 

As a simple model of the mirrors used in the experi- 
ments, let us consider a metallic slab with a large width, 
that is a width £ larger than a few plasma wavelengthes. 
We use expression (l§U|) of the force written as an integral 
over imaginary values of the frequencies (to — if, f real). 
Hence, the phase factor corresponding to one round trip 
inside the slab is a decreasing exponential with a real ex- 
ponent e~ 2Klt . For the plasma model (|14|> . n\ is given by 

ylf+ $+* and it is larger than £ for all values of 
f and k. When relaxation is taken into account, this is 
still the case except in a very narrow domain with values 
of f and k both close to zero. This domain has a negli- 
gible contribution to the integral H90J1 and it follows that 
the reflection amplitude of the slab may be replaced by 
the limiting expression obtained for the bulk. One thus 
recovers the Lifshitz expression for the Casimir force 0] 
which is widely used for comparing experimental results 
with theoretical expectations [l9j . 



D. Constraints on the force 

We now deduce general constraints which invalidate 
proposals made for tayloring the Casimir force at will by 
using specially designed mirrors ^3>EH- This generalizes 
to 3D space the results obtained for ID space in to 
which the reader is referred for further discussions. 

Expression (|90(l is an integral over the axis of imagi- 
nary frequencies essentially determined by the reflection 
amplitudes n [if] and ri [if] for f real. These amplitudes 
always have a modulus smaller than unity, for arbitrary 
dielectric multilayers (see appendix [UJ . They are nega- 
tive for arbitrary dielectric slabs (see appendix ^| and 
we deduce from the composition law l|47l) that this is still 
the case for arbitrary dielectric multilayers. It follows 
that the product of the reflection amplitudes of the two 
mirrors is always positive with a modulus smaller than 
unity 



o < n [if] r 2 [if] < 1 



(91) 



From this, we deduce first that the Casimir force has 
an absolute value smaller than the value (JTJ reached for 
perfect mirrors and that it remains attractive 



< F < F Ca 



(92) 



We also derive that the Casimir force decreases as a func- 
tion of the length 



OF 
dL 



< 



(93) 



This means that the properties obeyed by real mirrors 
strongly constrain the possible variation of the Casimir 
force, contrarily to what might have been expected at 
first sight mil. 

Note that we have considered mirrors used in the ex- 
periments which have electric permittivity but no mag- 
netic permeability. Different results would be obtained 
with magnetic mirrors, precisely with one of the mir- 
rors dominated by electric response and the other one 
by magnetic response. The product of the two reflection 
amplitudes would indeed be ne gati ve in this case and the 
Casimir force repulsive [6^, I63ll64l l65l | . 



VI. CONCLUSION 

We have presented a derivation of the Casimir force be- 
tween lossy mirrors characterized by arbitrary frequency 
dependent reflection amplitudes, in the Casimir geometry 
where the cavity is made with two parallel plane mirrors. 

We have shown how mirrors and cavities may be dealt 
with by using a quantum theory of optical networks. 
We have deduced the additional fluctuations accompa- 
nying dissipation from expressions of the optical theorem 
adapted to quantum network theory. The optical theo- 
rem is equivalent to the unitarity of the whole scattering 
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process which couples the modes of interest and the noise 
modes and it ensures that the quantum commutators of 
the output fields are the same as those of the input fields. 
The situation is different for the cavity fields which do not 
freely propagate. We have given a general proof of a the- 
orem previously demonstrated in particular cases |42l l4^ | 
which states that the modification of the commutators is 
determined by the usual Airy function, that is the spec- 
tral density associated with the Fabry-Perot cavity. For 
arbitrary lossy mirrors, the spectral density is determined 
by the reflection amplitudes as they are seen by the intra- 
cavity fields. It determines the radiation pressure exerted 
by vacuum fluctuations upon the mirrors with repulsive 
and attractive contributions associated respectively with 
resonant or antiresonant frequencies. The Casimir force 
is then obtained as an integral over the whole axis of 
real frequencies, including the contribution of evanescent 
waves besides that of ordinary waves. It is equivalently 
expressed as an integral over imaginary frequencies. The 
derivation only uses a few general assumptions certainly 
valid for real optical mirrors, namely causality, unitarity, 
high-frequency transparency for each mirror and stabil- 
ity of the compound cavity-vacuum system. It leads to 
a finite result without any further reference to a regular- 
ization technique 

The formula obtained in the present paper for the 
Casimir force was already known Du t, its scope of 
validity is widened by the present demonstration. It has 
been used to discuss the effect of imperfect reflection for 
the metallic mirrors used in the experiments. Different 
descriptions of the optical response of metals have been 
used, from the crude application of the plasma model 
(|14fl to a more complete characterization of the dielectric 
constant derived from tabulated optical data and dis- 
persion relations. This kind of calculations, discussed in 
great detail for the mirrors corresponding to the recent 
experiments (see for example [3lJ), has not been repro- 
duced here. Instead, we have presented general results 
valid for any real mirrors obeying the physical properties 
already evoked and shown that they strongly constrain 
the variation of the Casimir force. 

In the present paper, we have restricted our atten- 
tion on the limit of zero temperature although our work 
was partly motivated by a recent polemical discussion 
of the effect of temperature on the Casimir force be- 
tween real mirrors [llllllslllllllSllllllliilill- 
Since contradictory results may have raised doubts about 
the validity and consistency of various derivations of the 
Casimir force, we have considered it was important to 
come back to the first principles in this derivation. This 
has been done in the present paper for the case of zero 
temperature. A follow-on publication will show how to 
include the effect of thermal fluctuations in the treatment 
in order to obtain an expression free from ambiguities for 
the Casimir force between arbitrary lossy mirrors at non 
zero temperatures. 



APPENDIX A: THE DIELECTRIC SLAB 

In this appendix, we discuss in more detail the specific 
case of the dielectric slab. We consider lossy as well as 
lossless slabs. 

For a lossless dielectric medium, the permittivity e is 
real at real frequencies. For ordinary waves, kq and K\ 
are purely imaginary, so that the impedance ratios are 
real for both polarizations. Hence j3 is real (j3 — j3 r ) 
and a purely imaginary (a = ion) so that the scattering 
amplitudes l|44|) are read as 



t = 



sinh (3 r 



sinh (p r + ioti) 

sinh [3 r 



sinh [3 r cos on + i cosh p r sin on 
sinh (ioti) 



sinh (p r + ion) 

i sin a,- 



sinh (3 r cos cti + i cosh p r sin on 



(Al) 



The sum of the squared amplitudes is unity |i| 2 + |r| 2 = 1 
while the reflexion and transmission amplitudes are in 
quadrature to each other tr* +rt* = 0, which means that 
S is a unitary 2x2 matrix, as it was expected for a lossless 
mirror. This implies that the reflection amplitude has a 
modulus smaller than unity |r| < 1. This property also 
holds for lossy mirrors thanks to positivity of dissipation 
(see appendix Q. 

Unitarity is defined without ambiguity only in the case 
of ordinary waves. For a lossless slab and evanescent 
waves, kq is real - it is just the inverse of the penetration 
length of evanescent wave in vacuum - whereas n\ re- 
mains purely imaginary. Hence /3 as well as a are purely 
imaginary and it is no longer possible to obtain general 
bounds for the scattering amplitudes 



t 



sinh (i{3i) 



sin (3i 



sinh (i(3i + icti) sin (Pi + oti) 
sinh(iai) sin a.; 

sinh (ifti + icti) sin (Pi + on) 



(A2) 



In particular, |r| does not remain always smaller than 1 
(see more explicit discussions in appendix 151 with differ- 
ent results for the TE and TM polarizations). 

For imaginary frequencies finally, P and a are positive 
real numbers, for lossy as well as lossless slabs. In this 
case, general bounds are easily obtained for the ampli- 
tudes 



< t = 



sinh (P r 



< -r = 



sinh (p r + a r ) 
sinh (a r ) 



sinh (p r + a r ) 



< 1 
< 1 



(A3) 



The fact that r is negative with a modulus smaller than 
unity plays an important role in the derivation of con- 
straints on the Casimir force. 
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Interesting results are also obtained for the eigenvalues 
of the S— matrix, which have a simple form s± — r ± t 
since the slab is symmetrical in the exchange of its two 
ports. In the sector of ordinary waves, unitarity l|5U|l has 
a simple form in terms of s± = r ± i and of the similar 
quantities s'_j_ = r' ± t' defined on the noise matrix S' 

| s± | 2 + |4| 2 = l (A4) 

For the lossless slab, s± have a unit modulus and s± 
vanish. For a lossy slab, we have 

|s±| 2 < 1 (A5) 
This can be considered as a consequence of l|A4|) with 

l 1 2 

\ s ±\ ^ 0- Equivalently, it can be considered that uni- 
tarity i|A4jl fixes the modulus of s' ± when the modulus of 
a± is known. 

Condition i|A5(l will be found in appendix IU1 to express 
a passivity property for the slab, here for ordinary waves. 
This property still holds in the sector of imaginary fre- 
quencies, as a consequence of l|A3|l and of the following 
inequalities obeyed for all positive real numbers a and (3 

|sinh P =F sinh a\ < sinh (a -I- (3) (A6) 

Using the terms of appendix El this means that the do- 
main of passivity always includes the sectors of ordinary 
waves and imaginary frequencies, in the case of a dielec- 
tric slab. However, it does not necessarily include the 
sector of evanescent waves (see appendix |Bj . 

APPENDIX B: THE SECTOR OF EVANESCENT 
WAVES 

Ordinary waves correspond to frequencies u> > c |k| 
and real wavevectors k z whereas evanescent waves corre- 
spond to frequencies u> < c|k| and imaginary values of 
k z . Causal scattering amplitudes can be extended from 
ordinary to evanescent waves, by an analytical continua- 
tion through the physical domain of complex frequencies 
u> = i£ with > and 5Rk > 0. The 'energy condi- 
tions' which bear on quadratic forms are not necessarily 
preserved in this process. 

In order to illustrate the idea, let us consider the re- 
flection amplitude (|34|l at the interface between vacuum 
(eq = 1) and a lossless dielectric medium (ei real for u> 
real). In the sector of evanescent waves, k\ is imaginary 
and kq real, so that r and r are complex numbers with 
a unit modulus, that is also pure dephasings correspond- 
ing to the phenomenon of total reflection. Meanwhile, 
the transmission amplitudes differ from zero, which de- 
scribes how evanescent waves in vacuum are fed by the 
fields coming from the dielectric medium with an inci- 
dence angle larger than the limit angle. In these condi- 
tions, it is clear that the condition |r| 2 + |i| 2 < 1 fails. 

For the TE polarization, it turns out that 

|r TE | < 1 (Bl) 



in the evanescent sector at the interface between vac- 
uum and any dielectric medium. This property is always 
true in the sectors of ordinary waves and imaginary fre- 
quencies for an arbitrary mirror (see the appendices El 
and [UJ . Using high frequency transparency, it follows 
from the Phragmen-Lindelof theorem |66( that inequal- 
ity IjBlfl holds in the whole physical domain in the com- 
plex plane. This ensures that the closed loop function 
/ TE is analytic and, in particular, has no pole in the do- 
main > 0. In other words, since the open loop gain is 
smaller than unity, the closed loop cannot reach the os- 
cillation threshold, leading to the stability property used 
in the derivation of the Casimir force. 

Although it seems quite natural, this argument is not 
valid in the general case. For metallic mirrors for exam- 
ple, the condition |r| < 1 is violated in the evanescent 
sector for TM modes. The reflection amplitude is even 
known to reach large resonant values at the plasmon reso- 
nances 67]. Of course, this does not prevent the stability 
property to be fulfilled : the Fabry-Perot cavity is in this 
case a stable closed loop built on an open loop exceeding 
the unit modulus but with a phase such that the oscilla- 
tion threshold is not reached. 

We stress again that the stability property is necessary 
in the derivation of the Casimir force since it entails that 
the closed loop function is properly defined in the evanes- 
cent sector. When the more stringent property \r\ < 1 is 
also obeyed, it follows from expression l|86|l that the Airy 
function, which has been defined with the significance of 
a positive spectral density on ordinary waves, remains 
positive in the evanescent sector. When the property 
\r\ < 1 fails, the Airy function can no longer be thought 
of as a spectral density in the whole physical domain, 
but this does not invalidate the derivation of the Casimir 
force. 



APPENDIX C: THE DOMAIN OF PASSIVITY 

In this appendix, we discuss the related but not iden- 
tical properties corresponding to positivity of dissipation 
and passivity. 

We consider an arbitrary mirror, that is a recipro- 
cal network connecting two vacuum ports. For ordinary 
waves, we define the power dissipated by the mirror 

* = R nt eL n - er^r*) + (eM - eg^S 1 *) 

= ( ein l K n ) - ( e ° ut l l e ° ut ) ( C1 ) 

where we have introduced row vectors conjugated to the 
column vectors 

<e out | = |e out > + <e in | = \e in ) f (C2) 

This power is positive as a consequence of unitarity 

7T = {e ia \ I - Si S \e hl ) 

= (e in \S' <i S'\e m ) > (C3) 
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which corresponds to the positivity of the matrix I — S 

V|e) (e\I-S*S\e) > (C4) 

where |e) represents arbitrary input fields. Positivity can 
also be expressed in terms of the eigenvalues £ of S'S 

det (S^S-U) =0 £>0 (C5) 

These eigenvalues are always real and positivity of dissi- 
pation is equivalent to the fact that they are smaller than 
unity 

I < 1 (C6) 

Passivity is a property directly related to positivity of 
dissipation but defined more generally for complex fre- 
quencies in the physical domain. In order to discuss it, 
we extend the matrix S from the sector of ordinary waves 
through the analytical continuation already discussed. 
We extend S' similarly with the complex conjugation 
cautiously defined since it involves complex frequencies : 
conjugation corresponds to £ — > £* and k — ► k* and it 
preserves the physical domain -R£ > , SRk > 0; the 
derivations performed for an amplitude in the domain 
5R£ > , 3£ < are thus translated to similar deriva- 
tions for the conjugated amplitude in the quarter plane 
5R£ > , 3£ > 0. 

Then, the domain of passivity of S is defined by the 
domain of £ for which I— S is a positive matrix (ea lC4|) 
that is also for which the eigenvalues £ of S^S are smaller 
than unity (ea lC6(l . An important feature of this prop- 
erty is that it is stable under composition : when two net- 
works A and B are piled up as in Figure the quadratic 
forms appearing in (|C4fl simply add up so that passiv- 
ity of the network AB follows from passivity of the two 
networks A and B. This is a special case of a general 
theorem 0] which states that networks built up with 
passive elements are passive. 



Passivity means that the eigenvalues l—£ of the matrix 
/ — S are both positive, which is equivalent to the 
following inequalities 

Tr (I - S^S) > det (I - S f S) > (C7) 

It may be written in terms of the scattering amplitudes 

|r| 2 + |r| 2 + 2|i| 2 <2 

\rr-t 2 \ 2 >\r\ 2 + \r\ 2 + 2\t\ 2 -I (C8) 

Passivity implies that the scattering amplitudes have a 
modulus smaller than unity 

\r\ < 1 |f| < 1 \t\ < 1 (C9) 

Conversely, the latter conditions are necessary but not 
sufficient for passivity. 

For a mirror symmetrical in the exchange of its two 
ports, a slab for example, the passivity conditions take 
the simple form \r ± t\ < 1. The results of appendix lAl 
thus entail that the domain of passivity always includes 
the sectors of ordinary waves and imaginary frequencies, 
for arbitrary slabs fea lA5() . Using the stability of passiv- 
ity under composition, we deduce that this is also the case 
for arbitrary multilayers. It follows that the reflection 
amplitudes always have a modulus smaller than unity 
for imaginary frequencies (£ real) 

|rR]|<l (CIO) 
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